Astronomy & Astrophysics manuscript no. na 


February 2, 2008 


(DOI: will be inserted by hand later) 





Sodium abundances in nearby disk stars 



o 
o 



> 
in 
cn 
in 
in 
o 

o 

or 
6 



J.R. Shi 1 ' 2 , T. Gehren 2 , and G. Zhao 1 ' 2 

1 National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, P. R. China 

2 Institut fur Astronomie und Astrophysik der Universitat Miinchen, Scheinerstrasse 1, D-81679 Miinchen, Germany 

Received date / Accepted date 

Abstract. We present sodium abundances for a sample of nearby stars. All results have been derived from NLTE statistical 
equilibrium calculations. The influence of collisional interactions with electrons and hydrogen atoms is evaluated by comparison 
of the solar spectrum with very precise fits to the Nal line cores. The NLTE effects are more pronounced in metal-poor stars 
since the statistical equilibrium is dominated by collisions of which at least the electronic component is substantially reduced. 
The resulting influence on the determination of sodium abundances is in a direction opposite to that found previously for Mg 
and Al. The NLTE corrections are about —0.1 in thick-disk stars with [Fe/H] ~ —0.6. Our [Na/Fe] abundance ratios are about 
solar for thick- and thin-disk stars. The increase in [Na/Fe] as a function of [Fe/H] for metal-rich stars found by Edvardsson 
et al. (1993 1 is confirmed. Our results suggest that sodium yields increase with the metallicity, and quite large amounts of 
sodium may be produced by AGB stars. We find that [Na/Fe] ratios, together with either [Mg/Fe] ratio, kinematic data or stellar 
evolutionary ages, make possible the individual discrimination between thin- and thick-disk membership. 
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1. Introduction 

The 23 Na nucleus contains an extra neutron that makes its syn- 
thesis deviate from that of the a-nuclei such as 24 Mg. So the 
abundance ratios [Na/Fe] and/or [Na/Mg] in stars are a poten- 
tial source of information concerning the history of Galactic 
nucleosynthesis. 

It is expected that sodium is made during carbon and neon 
burning in massive stars, and it would therefore be ejected by 
Type II supernova (SN II). Type la supernova, on the other 
hand, produce little sodium (Tsujimoto et al. I1995i . However, 
23 Na can be synthesized as a primary or secondary element. In 
the first case 23 Na would be produced directly in the carbon 
burning process, and the production rate would be indepen- 
dent of the metal content of the proceeding stellar generation. 
In the second process, where 23 Na is produced by an excess 
of neutrons, the production rate would depend on that neutron 
excess which in turn would be determined by the pre-existing 
metal abundance. In the latter case 23 Na would be expected 
to display an underabundance compared with Fe or Mg in the 
most metal-poor stars. Although significant progress has been 
made in nucleosynthesis calculations, theoretical SN yields of 
sodium have not yet converged to consistency among mod- 
els from different authors (Woosley & Weaver 1995, Limongi 
et al. 



119981 Umeda et al. 200(1. Also some synthesis of 
sodium in intermediate-mass stars can occur in the hydrogen- 
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burning shell through the neon-sodium cycle (Denissenkov & 
Denissenkova 1990 Mowlavi fT999l . 

The observational picture for this element is complicated. 
The early studies report that sodium scales with iron at all 
metallicities (e.g. the review of Timmes et al. 1995| and 
Goswami & Prantzos 2000). However, new studies do not sup- 
port this classical picture. For metal-poor stars, Pilachowski et 
al. {1996 1 analyzed a sample of 60 halo subgiants, giants and 
horizontal branch stars in the interval — 3 < [Fe/H] < — 1, and 
found a small [Na/Fe] deficiency of —0.17 dex in the mean and 
that bright field halo giants do not show the excess of sodium 
found in their globular cluster counterparts. Baumuller et al. 
( 1998) performed a NLTE analysis of sodium and reported a 
clear trend of decreasing [Na/Fe] with decreasing metallicity. 
Similarly, Stephens ( 1999) found that [Na/Fe] decreases as one 
goes from [Fe/H]= —1 to [Fe/H]= —2, and thus shows the 
theoretically expected behaviour of odd-Z elements. Carretta 
et al. (2000) and Gratton et al. (2003) also found that sodium 
is overdeficient in stars with [Fe/H] < — 1; this result was con- 
firmed by very recent work of Gehren et al. ( 2004). 

For disk stars, Edvardsson et al. (1993) used high- 
resolution, low-noise spectra to determine sodium abundace 
in 189 F and G stars. Their LTE analysis suggested that the 
[Na/Fe] ratio is solar and varies very little with metallic- 
ity, which was confirmed by Chen et al. (2000). Based on 
5682/5688 and 6154/6160 lines, Prochaska et al. J2000t mea- 
sured sodium abundance in their thick-disk stars and found that 
all the stars show mildly enhanced [Na/Fe] and there is a mild 



trend with metallicity (see also Gehren et al. 2004). 
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Concerning the most metal-rich stars in the galactic disk, 
the study by Edvardsson et al. ([1993 1 raised a number of new 
questions. They found sodium relative to iron versus [Fe/H] 
shows an increase for [Fe/H] > 0.0 (see also Feltzing & 
Gustafsson 1998 their Fig. 8). Is the "upturn" real? Large star- 
to-star scatter was encountered. Could this scatter to be reduced 
by using better abundance criteria? 

It is known that [Fe/H] distributions of the halo, thick- and 
thin-disk stars reveal considerable overlap, so it can not al- 
low unequivocal classification for many stars without reference 
to other criteria such as age or kinematics. Fuhrmann ( 1998 
2002 1 suggested that the thick-disk stars represent the princi- 
pal population for the early Galactic evolution, while the stel- 
lar halo contributes by only minute amounts (see also Ryan & 
Smith 2003 1. Considering the general kinematics, there may be 
a significant overlap of all populations (Gehren et al. 2004 1. 
Gratton et al. jl996l 120001 |5003 ( have shown that thick-disk 
stars are more [Mg/Fe] enhanced than the thin-disk stars at 
the same [Fe/H] value. This trend has also been reported by 
Fuhrmann J1998I see also Prochaska et al. 2000 and Reddy 
et al. 2003 1. Finally, Mashonkina et al. (2003) have given evi- 
dence that [Eu/Mg] ratio allows the individual discrimination 
of halo stars from thick-disk stars, and similar behaviour is 
found for the [Al/Mg] ratio by Gehren et al. ( I2004> . 

The present work is based on a sample of nearby stars and 
aims at exploring their [Na/Fe] abundance ratios applying full 
spectrum synthesis based on level populations calculated from 
the statistical equilibrium equations. In Sect. 2 we present the 
observational techniques, while the atmospheric models and 
stellar parameters are discussed in Sect. 3. NLTE line forma- 
tion is discussed in Sect. 4. The discussion is presented in Sect. 
5, and the conclusions are found in the last section. 

2. Observations 

Our approach towards a representative abundance investigation 
is aimed at analyzing a volume-complete sample of the F, G 
and early K dwarfs in the disk population, north of declination 
6 = —15° and within 25pc of the Sun (see Fuhrmann 2004] for 
details). The spectra of our samples were obtained through the 
years 1995 to 2000 by Klaus Fuhrmann with the fiber-coupled 
Cassegrain echelle spectrograph FOCES (Pfeiffer et al. 1998 1 
mounted at the 2.2m telescope of the Calar Alto Observatory. 
Only few of them were exposed to a 1024 2 24/j CCD with 
A/ A A ~ 40000, and their wavelength ranges are limited to 
4000-7000 A, whereas the standard configuration was a 2048 2 
15/z CCD that covered 4000-9000 A with A/AA - 60000. All 
stars were observed at least twice with a signal-to-noise ratio 
of at least S/N - 150 up to S/N - 400 (see Fuhrmann 1998 
2000 for details). 

3. Atmospheric models and stellar parameters 

3.1. Model atmospheres 

Our analyses are all based on the same type of atmospheric 
model, irrespective of temperature, gravity or metal abundance. 
We use line-blanketed LTE model atmospheres, generated as 



discussed by Fuhrmann et al. (I1997> . The main characteristics 
are: 

- the iron opacity was calculated with the improved mete- 
oritic value log epo = 7.51 

- opacities for metal-poor stars with [Fe/H] < —0.6 were cal- 
culated using a-element (O, Ne, Mg, Si, S, Ar, Ca and Ti) 
abundances enhanced by 0.4 dex 

- the mixing-length parameter l/H p was adopted to be 0.5 
(see Fuhrmann et al. U993> 

3.2. Stellar parameters 

For most of the stars we adopt the stellar parameters deter- 
mined by Fuhrmann (1998 2000 1, where the effective tem- 
peratures are derived from the wings of the Balmer lines; sur- 
face gravities are taken from the strong line wings of the Mg lb 
triplet and compared with values obtained from Hipparcos par- 
allaxes. Iron abundances are based on Fe II lines, and the mi- 
croturbulence velocities are estimated by requesting that the 
iron abundance derived from Fe II lines should not depend on 
equivalent width. The uncertainties for the temperature, sur- 
face gravity, metal abundance and microturbulence velocities 
are generally assumed to be ±80 K, 0.1 dex, 0.2 km s _1 and 
0.07 dex respectively (see Fuhrmann 1998 and 2000 for de- 
tails). For Procyon the parameters were taken from Korn et al. 
(2003 1, where the temperatures are again from the wings of 
Balmer lines; surface gravities are based on the HIPPARCOS 
parallaxes, and iron abundances are taken from NLTE calcula- 
tions. 

Following Fuhrmann's work, we employ the macroturbu- 
lence parameters £ in the radial-tangential form and adopt val- 
ues as described in Gray ( 1984 1. The projected rotational ve- 
locity is obtained as the residual to the observed line profiles, 
after having removed the known instrumental profile obtained 
from the Moon spectra. 

4. NLTE line formation 

4. 1. Atomic model 

The model of Na we adopt goes back to the atomic proper- 
ties documented in Baumiiller et al. ([1998 1, where most of the 
problems have been described in detail. In this paper, we ex- 
tend the number of levels to n = 15. The model is closed with 
the Na + ground state. Fine structure splitting is neglected with 
the exception of 3p 2 Py 2 and 3p 2 P^ 2 . Part of that model was 
shown by Baumiiller et al. {1998 see their Fig. 2). Energy lev- 
els are taken from Martin & Zalubas (1981) if available. For 
levels with I > 4 hydrogenic approximations are used. All 
bound-free radiative cross-sections have been included from 
close-coupling calculations of Butler (1993 1 and Butler et al. 
(19931. 

To establish the influence of deviations from LTE on the 
ionization equilibrium of sodium, we have chosen the reference 
stars of Korn et al. (20031. The full analysis of the solar spec- 
trum (Kurucz et al. 1984 1 and those of the reference stars allows 
a reasonable choice of the hydrogen collision enhancement fac- 
tor resulting in Sh = 0.05 for Na , found by iteration together 



J.R. Shi et al.: Sodium abundances in nearby disk stars 



3 



Table 1. Atomic data of Na lines, log gf values and damping 
constants have been determined from solar spectrum fits 



A [A] 



Transition 



log 9f log c 6 



5682.642 
5688.214 
5889.959 
5895.932 
6154.228 
6160.751 
8183.260 
8194.800 



3p 2 P? /2 


-Ad 


2 D3/ 2 
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-30.05 


3p 2 P§ /2 


- 4d 


2 D 5/2 


-0.47 


-30.05 


3s 2 S 1/2 
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2p0 
r l/ 2 


+0.11 


-31.60 


3s 2 S 1/2 
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2 p0 
^3/2 
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-31.60 


3p 2 P? /2 


- 5s 


Sl/2 


-1.57 


-29.78 


3p 2 P§/ 2 


- 5s 


2 Si/ 2 


-1.28 


-29.78 


3p 2 P?/ 2 


- 3d 


D3/2 


0.28 


-30.73 


3p 2 P§/ 2 


- 3d 


2 D 5/2 


0.49 


-30.73 



with all basic stellar paraeters. This number is the same as in 
the previous work of Baumuller et al. ( 1998 1. 

All calculations have been carried out with a revised ver- 
sion of the DETAIL program (Butler & Giddings 1985 ) using 
accelerated lambda iteration (see Gehren et al. 2001 2004|for 
details). 

4.2. Atomic line data 

As our primary choice we take the oscillator strengths from 
the NIST data base. We note that the gf values are similar 
to those adopted by Takada et al. (2003), the average differ- 
ence log g/(this work)— log g/(Takeda) is —0.02 ± 0.04. For 
the line pair of 6154/6160, our gf values are same as the val- 
ues adopted by Chen et al. (2000), but slightly higher than the 
values of Edvardsson et al. ( 1993) and Feltzing & Gustafsson 
( 1998), the difference is about 0.04. However, in our present 
analysis the absolute value of the oscillator strengths is unim- 
portant because the abundances are evaluated in a fully differ- 
ential way with respect to the sun. 

Collisional broadening through van der Waals interaction 
with hydrogen atoms is important for strong Na D lines. As 
already pointed out by Gehren et al. (I2001ll2004t the resulting 
values of the van der Waals damping constants are mostly near 
those calculated according to Anstee & O'Mara's ( 119911119951 
tables. Fitting the solar profiles with a common Na abundance 
then shows that our damping constants for the resonance lines 
are roughly a factor of 2 below those calculated. The adopted 
line data {gf values and van der Waals damping constants) are 
given in Table 1 . 

4.3. Sodium abundances 

The abundance determinations for our program stars are made 
using spectral synthesis. The synthetic spectra are convolved 
with macroturbulence, rotational and instrumental broadening 
profiles, in order to match the observed spectral lines. Based 
on 5682/5688, 5890/5896 and 6154/6160 line pairs, our abun- 
dance results are obtained from profile fits. However, for Na D 
lines, as pointed out by Baumuller et al. {1998 ), even in the so- 
lar atmosphere the inner cores of the lines are affected by devi- 
ations from LTE. We note that even for some metal-poor stars, 
only the line wings can be fitted in LTE. So the sodium abun- 
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Fig. 1. Comparison of derived [Na/Fe] for stars in common 
with other studies. Filled circles (•) are from Baumuller et al. 

; as- 



( 1998 ); crosses ( x) are from Feltzing & Gustafsson ( 1998 ) 
terisks (*) are from Carretta et al. (2000); pluses (+) are from 
Fulbright (2000); diamonds (o) are from Gratton et al. (2003 ); 
open circles (o) are from Takeda et al. ( 2003 ); triangles (A) are 
from Reddy et al. d2003l 



dances are obtained until a best fit to the observed line wings 
is provided for the LTE analyses. Considering the NLTE abun- 
dances, out results do not show large abundance differences 
between different lines (however, see Takeda et al. 2003 ). The 
final abundance scatter of single line is between 0.01 and 0.06 
with < cr[Na/Fe]>= 0.020 ± 0.007 for all stars. The derived 
abundances are presented in Table 2. Complete results for in- 
dividual lines of all stars can be obtained from the author on 
request. 

4.4. Comparison with other works 

Sodium abundances have been determined by several groups 
including both LTE and NLTE analyses. In Fig.l we compare 
the [Na/Fe] values (NLTE) determined in this paper with those 
from the literature. The agreement is generally good. However, 
a closer inspection reveals some systematic differences. In the 
remaining part of this section we will briefly discuss some pos- 
sible reasons for these differences. 



Baumuller et al. (1998 



Baumuller et al. ( 1998 1 performed a full NLTE line for- 
mation for a sample of metal-poor stars. We have four stars in 
common with theirs. Our NLTE results are very much in agree- 
ment with theirs, as shown in Fig. 1. The average differences 
between our [Na/Fe] and theirs are 0.02 ± 0.10. We note that 
the considerable scatter is mostly due to the different stellar 
parameters adopted. 

Feltzing & Gustafsson ( 1998) 

In this study of 47 G and K metal-rich stars, they confirmed 
the increase of [Na/Fe] as a function of [Fe/H] found previously 



by Edvardsson et al. (1993). Their results were determined 
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Table 2. Atmospheric parameters and sodium abundances of Table 2 (continued) 
the program stars. Rotational velocities and macroturbulence 

data are the same as in Fuhrmann ( 1998l l2000t ™> T ^ lo S5 [ Fe/H J £ [Natfe] [Na/Fe] 
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from 6154/6160 lines, the gf values adopted in this study are 
slightly lower than ours, the differences are —0.02/ — 0.05, re- 
spectively. The sodium abundance determined by this study is 
in good agreement with ours. For three stars in common, the 
average difference is 0.01 ± 0.04. 

Carretta et al. (2000> 

This work reanalysed the data of Edvardsson et al. ( 1993 1 
including a NLTE correction. They found no large differences 
from the LTE results of Edvardsson et al., which is shown in 
their Fig. 9. The gf values of the 6154/6160 lines adopted 
by Edvardsson et al. dl993> are slightly lower than ours, the 
differences are —0.04 and —0.03, respectively. The results of 
Carretta et al. (2000 1 are very much in agreement with ours. 
For the 46 common stars from Edvardsson et al. (Q993 1, the 
average difference between ours and theirs is —0.01 ± 0.05. 

Fulbright(2000 2002 1 

This analysis deals with a large number of metal-poor 
stars, of which five stars are in common with our sample. 
As pointed out by Gehren et al. (2004i, there are large dif- 
ferences in stellar parameters; on average, the typical differ- 
ence in temperatures for the individual star is about 150 K, 
and there are a few stars for which differences may be much 
larger. Fulbright's abundances are slightly larger than ours (see 
Fig. 1). For the five stars in common with our list we ob- 
tain A[Na/Fe]=-0.09 ± 0.05. We find that the largest differ- 
ence (^0.4 dex) comes from warm metal-poor stars such as 
HD 284248. We note that considerable discrepancies remain 
even after correcting for the difference in effective tempera- 
tures. As these stars are metal-poor, they may suffer from rela- 
tively large NLTE effects (see Sect 5.1). 

Gratton etal.(2003) 

The authors considered NLTE corrections and reported the 
sodium abundances for 150 field subdwarfs and subgiants. Our 
results are mostly in agreement with theirs. For the six stars in 
common, the average difference is 0.07 ± 0.05. 

Takeda et al. (2003 1 

Very recently, Takeda et al. ( 2003 1 performed a NLTE anal- 
ysis of sodium abundances for disk and halo stars. They re- 
analysed the data from Chen et al. ((2000 1. Their results vary 
from line to line (see their Tables 4 and 5). The results for 
the 6154/6160 lines show no large differences from those of 
Chen et al. (see their Fig. 5c). Our [Na/Fe] values are in agree- 
ment with theirs for the lines of 6154/6160 A for the three stars 
in common; the average difference is 0.04 ± 0.06. As already 
pointed out by Takeda et al., this line pair is not very sensitive 
to the NLTE effects. 

Reddy et al. (2003) 

This sample consists of 187 F and G dwarfs. The sodium 
abundances were based on 6154/6160 lines, the gf values 
adopted in their work are nearly same as those we adopted. 
A very close star by star agreement between the [Na/Fe] ratios 
derived here and those measured by Reddy et al., no differ- 
ence exceeds 0.20 dex. For the six stars in common, we obtain 
A[Na/Fe]=0.00±0.08. The large scatter is due to HD 218470 



(0.16 dex). This is a high vs'mi (13.1 km s _1 ), high tempera- 
ture (T c ff = 6407) and low gravity (log g = 4.02) star. For this 
star their temperature is about 70 K higher than ours. 

5. Discussion 

5. 1. NLTE effects 

The abundance analyses of sodium clearly show the NLTE ef- 
fect. There is a tendency that the NLTE effect is large for warm 
metal-poor subgiant stars, as would be expected. Baumiiller et 
al. 0998 1 show in their Fig. 11 that the NLTE line cores for 
the Na D lines are much deeper than in LTE in metal-deficient 
stars. In LTE abundance analyses this is compensated by sim- 
ply increasing the sodium abundance until the observed equiv- 
alent width is reproduced. Such results are displayed in Table 3 
of Baumiiller et al. (see also their Table 2). The results confirm, 
starting from the solar abundance towards the lower metallicity, 
that the NLTE abundance effect resulting from the Na D lines 
increases reaching a maximum near [Fe/H]= —2. LTE abun- 
dances can be significantly different from their NLTE counter- 
parts, with differences reaching more than 0.6 dex in extreme 
cases. In Table 3 it is evident that the NLTE effects are system- 
atically stronger for the hotter models, which is in agreement 
with the statistical equilibrium of both aluminium and magne- 
sium (Baumiiller et al. 119971 Zhao & Gehren |2000T >. As ex- 
pected, the strongest departures from LTE are found for mod- 
els with high temperature and low metallicity. The reduction 
of surface gravity results in a decreased efficiency of colli- 
sions with electrons and hydrogen atoms, which again leads 
to stronger NLTE effects. In Fig. 2, the differences between 
LTE and NLTE analyses for our program stars are plotted as a 
function of metal abundance, temperature and luminosity. 

It should be noted that the NLTE effects differ from line 
to line, reflecting their individual properties. For example, the 
8183/8194 lines show relatively large NLTE effects compared 
to the other lines at similar line-strength. From Table 3, we 
can see the importance of the NLTE effect on different lines 
in sodium abundance determination: 

- the weak 6154/6160 lines show the smallest NLTE abun- 
dance effects (< 0.1 dex) for the disk stars, so they are the 
best abundance indicators for LTE analyses of moderately 
metal-poor stars. 

- the NLTE corrections for 5682/5688 lines are relatively 
small (~ 0.1 dex) for dwarfs; but they increase for sub- 
giants (> 0.1 dex). Therefore these lines should no longer 
be analyzed with LTE. 

- for the 5890/5896 lines, which are important for determin- 
ing the sodium abundances of very metal-deficient stars, 
it is necessary to take account of the NLTE effects except 
for very cool and metal-rich dwarfs. Especially for warm 
metal-poor giants, the NLTE correction can reach 0.6 dex 
(see Table 3 for details). 

- the 8183/8194 lines are also important for determining the 
sodium abundances of very metal-deficient stars because 
of their strength. We suggest that, even for cool metal-rich 
stars, the NLTE effects should be taken into account. 
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Fig. 2. Difference of [Na/Fe] abundance ratios calculated under 
NLTE and LTE assumptions as a function of metal abundance 
(a), temperature (b) and luminosity (c). Filled circles refer to 
the thin-disk stars, open circles refer thick-disk stars. 



A warning is appropriate that Table 3 should not be used 
for LTE abundance corrections if precise Na abundances are 
aimed at. The data may vary considerably whenever the mi- 
croturbulence velocity deviates from 1 kms -1 , and any simple 
correction may result in an increased scatter. The reference to 
[Na/Fe] = may also be a bad approximation, particularly in 
metal-poor stars. 

5.2. Sodium abundance and nucleosynthesis in the 
Galaxy 

The variation of [Na/Fe] with the stellar metal abundance 
[Fe/H] contains information about the chemical evolution of 
the Galaxy. Fig. 3 displays the run of [Na/Fe] ratio (calculated 
in NLTE) with the overall metal abundance for all stars con- 
sidered in this paper. One important feature that can be seen 
from Fig. 3 is that the bulk of the sodium abundances scale 
approximately as iron in thin-disk stars, with a small overabun- 
dance for the thick-disk stars. The "upturn" in [Na/Fe] versus 
[Fe/H] for metal-rich disk stars, as observed by Edvardsson et 
al. £l993l, is also reproduced 

Using magnesium instead of iron as the reference can re- 
move Type la SNe from consideration, because nearly all 
the sodium and magnesium are produced in massive stars 
(Timmes et al. 11995b . Calculations of nucleosynthsis for SNe 
la show that they should produce only very little amounts of 
sodium (Nomoto et al. 11997b . The odd-even effects will be 
enhanced using magnesium as the basis of comparison. The 
overall behaviour of [Na/Mg] ratios versus [Mg/H] is there- 
fore shown in Fig. 4 where the Mg abundances are taken from 
Fuhrmann fl998l 12000). We confirm Arnett's ( 11971b predic- 
tion that [Na/Mg] ratios decline with decreasing metallicity. 
This suggests that sodium abundances may be more closely 
coupled to chemical evolution history, with a gradual enrich- 
ment of [Na/Mg] from thick to thin disk. 

Presently, the source of sodium remains to be confirmed. 
Thus it is difficult to use sodium as a probe of Galactic chemi- 
cal evolution until the effects of individual stars can be quanti- 
fied. Theoretical sodium yields of SN II have not yet converged 
to consistency among models of different authors (Woosley & 
Weaver [19951 Umeda et al. 120001 Umeda et al. ( 120001 show 
the yields of sodium are much less for Z = compared to 
the Z = 0.02 for a 20 Mq model, while Woosley & Weaver 
(11995b predict that for a 22 Mq SN model the amount of 
sodium decreases with increasing metallicity in the range of 
Z = to Zq. Also, AGB stars are potential sites for the pro- 
duction of primary sodium (Mowlavi 11999b . In this scenario 
synthesis of sodium can occur in hydrogen-burning shells dur- 
ing ON-cycling (Denissenkov & Tout I2000T > . This process is 
supported by observations of red giant stars in globular clus- 



ters that appear to be sodium-rich (Kraft et al. 1997 Gratton et 
al. l200TV 



We confirm that departures from LTE of the Na level popu- 
lations appear to be larger for warm and giant stars, which was 
also found by Baumiiller et al. ( 1998 1 and Takeda et al. (2003 1. 



1 We have 47 stars common with the work of Edvardsson et al. 
(19931 . We note that more metal-rich stars are included in our sam- 
ple. The upturn in [Na/Fe] vs. [Fe/H], as observed by Edvardsson et 
al. is reproduced and reinforced. We confirm that the upturn is real 
(also see Feltzing & Gustafsson 1998 1. 
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Table 3. Abundance differences between NLTE and LTE obtained by fitting calculated LTE equivalent widths of Na I lines using 
the same parameters for various stellar model atmospheres based on the line-blanketed model grid. Results refer to A[Na/Fe] = 
[Na/Fe]NLTE~ [Na/Fe]LTE- No entries are given for extremely weak lines. The calculations all refer to £ = 1 km s _1 and [Na/Fe] 
= 



r efi logg [Fe/H] 5682 5688 5890 5896 6154 6160 8183 8194 A [Na/Fe] 
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Fig. 3. Abundance ratios [Na/Fe] for NLTE analyses. The Fig. 4. [Na/Mg] abundances ratios for NLTE analysis as a func- 



meaning of the symbols is the same as in Fig. 2 



tion of [Mg/H]. Symbols are the same as in Fig. 2 



Detailed modelling of the Galactic chemical evolution has 
been attempted by many authors (e.g. Timmes et al. 119951 
Goswami & Prantzos 2000 1. Based on Woosley & Weaver's 
dl995> metallicity-dependent yields, the first authors calculate 
the behaviour of [Na/Fe] as a function of metallicity. Their 
results predict that [Na/Fe] decreases from [Fe/H] ~ —3 to 



~ —1.5, while it increases from [Fe/H] ~ —1.5 to ~ —0.5 (see 
their Fig. 17). However using the same yields with an iron yield 
reduced by a factor of 2, Goswami & Prantzos' calculation ex- 
plains the decline of [Na/Fe] from [Fe/H] 1 to [Fe/H] 2 

reasonably well. Their prediction suggests an ever-decreasing 
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[Na/Fe] toward the lower metallicity regime. They used a dif- 
ferent initial mass function (IMF), and a different halo model. 

Our observational results provide some implications for the 
nuclesynthesis of sodium. We note that in halo stars sodium 
is significantly underabundant relative to iron, and [Na/Fe] (or 
[Na/Mg]) increases with metallicity (Fig. 5). In the region of 
overlapping metal abundances [Na/Fe] shows a clear distinc- 
tion between halo and thick-disk stars. The linear relationship 
indicates that sodium production is connected to the metallicity 
of the progenitor SN as suggested by Nissen & Schuster ( 1997 1. 
This result agrees with nucleosynthesis calculations of massive 
stars by Umeda et al. yOOO), but it is at variance with the cal- 
culations by Woosley & Weaver ( 1995 1. We note that there are 
no extremely metal-poor stars in the sample of Gehren et al. 



(20041, whereas Mc William et al. (11995 ) find that, for stars 
with [Fe/H] < —3, [Na/Fe] rises with decreasing metal abun- 
dance. Based on this result Tsujimoto et al. (2002i argued that 
sodium yields decrease with increasing Z for the very metal- 
poor stars. However, nearly all the stars are giants and the abun- 
dances were derived from the Na D (5890/5896) lines, which 
are very sensitive to NLTE effects (see Sect. 5.1; Takeda et al. 
2003 1. There are two possibilities to explain the high [Na/Fe] 
values of Mc William et al., 



- NLTE corrections for Na D lines are large for metal-poor 
cool giants, so if it is considered, the [Na/Fe] ratios would 
show an abundance deficiency (however, see Takeda et al. 
120031 . 

- sodium is also produced in metal-poor giant stars (Goriely 
& Siess 2001 ), which could explain why the Pb-stars also 
have a large sodium content (Van Eck et al. l2001i . 



A definitive confirmation of the latter assumption will require 
more detailed nucleosynthesis calculations. It is highly desir- 
able to include the sample of ultra-low metallicity stars for es- 
tablishing the behaviour of [Na/Fe] in this region, which is im- 
portant for investigating the history of the early Galaxy as well 
as for constructing the model of Galactic chemical evolution. 

Based on [Eu/Ba] ratios, Mashonkina et al. (2003 1 sug- 
gested that during the active phase of the thick disk forma- 
tion evolved low mass stars may enrich the interstellar gas. 
Our observational results show that [Na/Fe] values are about 
solar, possibly being slightly overabundant in the thick-disk 
stars, which would confirm the suggestion that AGB stars 
may contribute some additional sodium. However, metallicity- 
dependent SN yields of sodium can not be excluded. 

While the origin of the "upturn" for the metal-rich stars 
is still unclear, Feltzing & Gustafsson (1998 1 argued that 
metallicity-dependent SN yields may play a role. Though a 
rather insufficient fraction of stars supersolar metal abundances 
prevents us from making any definitive argument, we sug- 
gest that quite large amounts of Na may be produced by pro- 
ton capture on 22 Ne nuclei for AGB stars (Denissenkov & 
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Fig. 5. Top: abundance ratios [Na/Fe] for NLTE analysis as 
a function of [Fe/H]. Bottom: abundance ratios [Na/Mg] for 
NLTE analysis as a function of [Mg/H]. Filled and open circles 
refer to the thin- and thick-disk stars in this paper, respectively. 
Asterisks refer to the stars from Gehren et al. 20041 



5.3. Abundance ratios, kinematics, ages and the 
Galactic stellar populations 

It is important to distinguish the individual membership of 
stars in our sample. We are aware that such a discussion re- 
quires more than the interpretation of abundance ratios alone. 



Following Fuhrmann's (1998 2000) work we use [Mg/Fe] 
= 0.25 as a first step to discriminate between thin- and thick- 
disk stars. 

5.3.1 . Abundance ratios and population 
membership 

Fuhrmann ( 1998 1 attributed the scatter in [Mg/Fe] to a mixing 
of thin- with thick-disk stars and a different chemical evolution 



Denissenkova ll99"0l . 



of these two stellar populations (see also Gratton et al. 2000 
2003 for remarks on O and Mg with respect to Fe). He sug- 
gested that the stars with higher [Mg/Fe] at a given [Fe/H] are 
old thick-disk stars (see also Reddy et al. l2003i . 
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Fig. 6. [Na/Fe] abundance ratios for NLTE analyses as a func- 
tion of [Mg/Fe]. The meaning of symbols is the same as in Fig. 
2 

We show the distribution of [Na/Fe] versus [Mg/Fe] for the 
two populations of thick- and thin-disk stars in Fig. 6. Even 
comparing abundance alone, Fig. 6 suggests differences among 
the two populations: the thin-disk stars are well separated on 
the left side. 

5.3.2. Kinematic properties and population 
membership 

The original classification of stellar populations was based 
much more on kinematics than on anything else. It is therefore 
of high priority to check if kinematic properties and abundance 
ratios have more in common than a coarse change with the 
mean metal abundance itself. Fuhrmann (2000 2002; see also 
Ryan & Smith 2003 1 argues that in fact halo stars may consti- 
tute an extremely small minority of all Galactic stars, whereas 
most of the stars with intermediate kinematics instead belong 
to the thick-disk. Fuhrmann (2000 1 also suggested that there 
is some kinematical overlap of the disk populations. Therefore 
the determination of Galactic kinematics provides an important 
test. 

Fuhrmann (2000 1 has discussed in more detail the kine- 
matical aspects of our sample. The basic data are available 
in a number of catalogues which are electronically accessible 
(e.g. via http://www.ari.uni-heidelberg.de/aricnsi. The correc- 
tion for the basic solar motion is that of Dehnen & Binney 
(Q998), namely, U /V /W = 10.00/5.25/7.17 km s" 1 . 
Whenever these data are not provided, we have calculated the 
space velocities following the model described by Chen et al. 
(120001 . 

The first inspection of Fig. 7 shows that the V velocities 
of thin-disk stars tend to be significantly larger than those of 
the other populations, although there is some overlap between 
thin- and thick-disk stars. The two most exceptional cases are 
the thin-disk star HD 5271 1 (V=-98 km s _1 ) and the transition 
star HD 90508 (V=-89 km s _1 ), and we will come back to 
them later. The correlation between low [Na/Mg] ratios and low 
V or high peculiar space velocities (v spec ) is also instructive 
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Fig. 7. Correlation between Na abundance ratios and orbital V 
velocity components. Symbols are as in Fig. 2 



(Fig. 8). The kinematic status in the Toomre diagram of our 
sample was shown by Fuhrmann (2000 in his Fig. 17). 

We now come back to the transition stars HD 90508 
(U/V/W = 31/ - 89/30 km s" 1 ) and thin-disk object 
HD52 711 (U/V/W = -17/ - 98/ - 13 km s" 1 ): Fuhrmann 
(20001 argues that although HD 90508 shows thick-disk kine- 
matics, its chemistry and age are more like thin-disk stars. 
HD 52711 is thought to be a thick-disk star based on its kine- 
matics (e.g. Ibukiyama & Arimoto 2002 1, however, it is very 
young (6.8 Gyr) and [Mg/Fe] is about solar (0.04), therefore it 
is also more like a thin-disk star. 

5.3.3. Stellar ages and population membership 

We now turn to the age aspects of analyzed stars. Gratton et al. 
( 2000 1 concluded that there is a gap in the distribution of [O/Fe] 
as expected from a sudden decrease in star formation during the 
transition from the thick- to thin-disk phase. Fuhrmann (2000 1 
suggested that the key feature to identify thick- and thin-disk 
stars individually is provided by the stellar ages, while there is 
a considerable age overlap between halo and thick-disk stars. 
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Fig. 8. Correlation between [Na/Fe] abundance ratios and pe- 
culiar space velocities v spec = (U 2 + V 2 + W 2 ) 1 ! 2 (top), 
and correlation between [Na/Mg] abundance ratios and pecu- 
liar space velocities v spe c — (U 2 + V 2 + W 2 ) 1 ! 2 (bottom). 
The meaning of the symbols is the same as in Fig. 2 



Using the stellar effective temperatures together with the 
absolute magnitudes based on Hipparcos parallaxes, masses 
and approximate ages can be interpolated according to [Fe/H] 
and [a/Fe] from adequate tracks of stellar evolution. Such cal- 
culations have recently become available through the work of 
VandenBerg et al. ( 120001 12002i and Yi et al. ( l2"5o"3l. In th is pa- 
per we adopt the Yonsei-Yale isochrones (Yi et al. 2003 1 since 
the evolutionary tracks of VandenBerg contain only metal-poor 
([Fe/H] < -0.31) and low mass stars (M stm < 1.20). The 
Yonsei-Yale isochrones are calculated with new OPAL opac- 
ities and Kurucz model atmospheres for a set of metallicities 
Z = 0.00001, 0.0001, 0.0004, 0.001, 0.004, 0.007, 0.01, 0.02, 
0.04, 0.06, 0.08, and [a/Fe] = 0.0, 0.3, 0.6. The mixing length 
and helium enrichment rate were fixed to a = I / H p = 1 . 743 1 
and AY/AZ = 2.0, respectively. The full set of stellar models 
and a FORTRAN package that works for mass, metallicity and 
a-enhancement interpolation are available from the authors. 

The ages derived here are based on the assumption that 
[0/Mg]= 0. From Fig. 9 we can see that all the thin-disk 



stars do not exceed ~ 9 Gyr, whereas with the exclusion 
of HD 165401 all thick-disk stars are older than — 9 Gyr. 
Fuhrmann (1998 1 identified this star as a thick-disk member 
with an age of about 10 Gyr, but in his paper II (Fuhrmann 
2000) he noted that it is a chromospherically active, high- 
velocity star, so it was excluded from his later analysis. 

As already pointed out by VandenBerg et al. ( 2002 1 and it- 
erated by Gehren et al. (2004), a significantly higher [Mg/Fe] 
ratio could reduce the ages of thick-disk and halo stars by as 
much as 3 Gyr. For comparison, we also calculated the ages 
for these stars without a-enhancement. Then our result is fairly 
in agreement with Fuhrmann's, namely that all the thick-disk 
(excluding HD 165401) and halo stars are older than 12 Gyr, 
while all the thin-disk stars are younger than 9.5 Gyr except the 
three transition stars (HD 90508, HD 143761 and HD 187923). 
Therefore Fuhrmann's result that the key feature to identify 
thick- and thin-disk stars individually is provided by the stel- 
lar ages still holds. However, the 3 Gyr long star formation gap 
between thin- and thick-disk is not confirmed here. 

The three transition stars (HD 90508, HD 143761 and 
HD 187923) have ages about 8.8 Gyr, thus they are as old as 
the oldest thin-disk stars. They would be about 10.5 Gyr old if 
no a-enhancement were considered. 

We also note that the age of the thick-disk star HD 6582 
seems to be very high (>15 Gyr); Gehren et al. (|2004) have 
discussed this in detail. 

6. Conclusions 

We have determined sodium abundances for some 90 nearby 
stars, spanning the range -0.9 < [Fe/H] < +0.4. All abun- 
dances are derived from NLTE statistical equilibrium calcula- 
tions. Based on our results we come to the following conclu- 
sions: 

1 . The [Na/Fe] ratios are about solar for thin- and thick-disk 
stars. The "upturn" of [Na/Fe] for metal-rich stars is repro- 
duced. 

2. The NLTE effects are different from line to line. The 
6154/6160 lines are comparatively less sensitive to NLTE 
effects, while the Na D and near-IR line pairs show large 
NLTE effects. Large departures from LTE appear for warm 
metal-poor subgiant stars. 

3. Our results suggest that the sodium yields increase with 
metallicity, and quite large amounts of sodium may be pro- 
duced by AGB stars. 

4. Combined with the abundance ratios, kinematics and stel- 
lar ages, we identify individual thick- and thin-disk stars. 
Fuhrmann's suggestion that stellar ages provide a key fea- 
ture in discriminating the two disk populations is supported. 

It would be important to perform NLTE abundance deter- 
minations for both Mg and Al for these program stars to con- 
firm our identification, because both Al and Mg show large 
NLTE effects for metal-poor stars (Gehren et al. 12004b . It will 
be also necessary to investigate some extremely metal-poor 
stars. 
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